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Abstract. A detailed tight-binding analysis of the electron band structure of 
the Cu02 plane of layered cuprates is performed within a a-band Hamiltonian 
including four orbitals - Cu3(ia;2_^2, Cu45, 02^^^, and 02py. Both the experimental 
and theoretical hints in favor of Fermi level located in a Cu or O band, respectively, 
are considered. For these two alternatives analytical expressions are obtained for the 
LCAO electron wave functions suitable for the treatment of electron superexchange. 
Simple formulae for the Fermi surface and electron dispersions are derived by 
applying the Lowdin down-fold procedure to set up the effective copper and oxygen 
Hamiltonians. They are used to fit the experimental ARUPS Fermi surface of 
Pbo.42Bii.73Sri.94Cai.3Cui.9208+cc and both the ARPES and LDA Fermi surface of 
'Nd2-xCexC\i04-s . The value of presenting the hopping amplitudes as surface integrals 
of ab initio atomic wave functions is demonstrated as well. The same approach 
is applied to the Ru02 plane of the ruthenate Sr2Ru04. The LCAO Hamiltonians 
including the three in-plane 7r-orbitals KuMxy^ Og,2py^ O^^Px and the four transversal 
TT-orbitals KuMzx^ KuMyz^ Oq2pz^ 0\^2pz^ are separately considered. It is shown that 
the equation for the constant energy curves and the Fermi contours has the same 
canonical form as the one for the layered cuprates. 
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1. Introduction 

After the discovery of the high- Tc superconductors the layered cuprates became one 
of the most studied materials in the solid state physics. A vast range of compounds 
were synthesised and their properties comprehensively investigated. The electron band 
structure is of particular importance for understanding the nature of superconductivity 
in this type of perovskites [|l[ . Along this line one can single out the significant success 
achieved in the attempts to reconcile the photoelectron spectroscopy data |^ and the 
band structure calculations of the Fermi surface (FS) especially for compounds with 
simple structure such as Nd2-xCea;Cu04_(5 P, |[. A qualitative understanding, at least 
for the self-consistent electron picture, has been achieved and for the most electron 
processes in the layered perovskites one can employ adequate lattice models. 

There is not much analysis of the electronic band structures of the high-Tc materials 
in the terms of single analytical expressions available. This is something for which there 
is a clear need, in particular to help in the construction of more realistic many-body 
Hamiltonians. The aim of this paper is to analyse the common features in the electron 
band structure of the layered perovskites within the tight-binding (TB) method 0. 
In the following we shall focus on the metallic (being eventually superconducting) 
phase only, with the reservation that the antiferromagnetic correlations, especially 
in the dielectric phase, could substantially change the electron dispersions. It is 
shown that the linear combination of atomic orbitals (LCAO) approximation can be 
considered as an adequate tool for analysing energy bands. Within the latter exact 
analytic results are obtained for the constant energy contours (CEC). These expressions 
are used to fit the FS of Nd2-xCexCvL04-s [i], Pbo.42Bii.73Sri.94Cai.3Cui.9208+x @], 
and Sr2Ru04 measured in angle-resolved photoemission/angle-resolved ultraviolet 
spectroscopy (ARPES/ARUPS) experiments. 

In particular, by applying the Lowdin perturbative technique for the Cu02 plane we 
give the LCAO wave function of the states near the Fermi energy . These states could 
be useful in constructing the pairing theory for Cu02 plane. For the layered cuprates we 
find an alternative concerning the Fermi level location — Cu3rfa^2_^2 vs. 02pa character 
of the conduction band. It is shown that analysis of extra spectroscopic data is needed 
in order for this dilemma to be resolved. As regards the Ru02 plane, the existence of 
three pockets of the FS unambiguously reveals the RnMs character of the conduction 
bands g|. 

To address the conduction bands in the layered perovskites we start from a common 
Hamiltonian including the basis of valence states 02p and RuAds^ or Cu3rfa;2_^2, CnAs 
respectively for cuprates. Despite the equivalent crystal structure of Sr2Ru04 |p!0|] 
and La2-a:Baa;Cu04 [|Tl|, the states in their conduction band(s) are, in some sense. 



complimentary. In other words, for the Cu02 plane the conduction band is of a-character 
while for the Ru02 plane the conduction bands are determined by tt valence bonds. This 
is due to the separation into a- and vr-part of the Hamiltonian H = + H^^"^ in first 
approximation. The latter two Hamiltonians are studied separately. 
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Accordingly, the paper is structured as follows. In section]^ we consider the generic 
Hamiltonian of the CuOs plane ig, [g and H^^^ = H^'^y^ + H^'^ is then studied 
in section The results of the comparison with the experimental data are summarised 
in section ^. Before embarking on a detailed analysis, however, we give an account of 
some clarifying issues concerning the applicability of the TB model and the band theory 
in general. 



1.1. Appology to the band theory 

It is well-known that the electron band theory is a self-consistent treatment of the 
electron motion in the crystal lattice. Even the classical 3-body problem demonstrates 
strongly correlated solutions, so it is a priori unknown whether the self-consistent 
approximation is applicable when describing the electronic structure of every new 
crystal. However, the one-particle band picture is an indispensable stage in the 
complex study of materials. It is the analysis of experimental data using a conceptually 
clear band theory that reveals nontrivial effects: how strong the strongly correlated 
electronic effects are, whether it is possible to take into account the influence of some 
interaction-induced order parameter back into the electronic structure etc. Therefore 
the comparison of the experiment with the band calculations is not an attempt, as 
sometimes thought, to hide the relevant issues — it is a tool to reveal interesting and 
nontrivial properties of the electronic structure. 

Many electron band calculations have been performed for the layered perovskites 
and results were compared to data due to ARPES experiments. The shape of the Fermi 
surface is probably the simplest test to check whether we are on the right track or some 
conceptually new theory should be used from the very beginning. 

The tight-binding interpolation of the electronic structure is often used for fitting 
the experimental data. This is because the accuracy of that approximation is often 
higher than the uncertainties in the experiment. Moreover, the tight-binding method 
gives simple formulae which could be of use for experimentalists to see how far they can 
get with such a simple minded approach. The tight-binding parameters, however, have 
in a sense "their own life" independent of the ab initio calculations. These parameters 
can be fitted directly to the experiment even when, by some reasons, the electron band 
calculations could give wrong predictions. In this sense the tight-binding parameters 
are the appropriate intermediary between the theory and experiment. As for the theory, 
establishing of reliable one-particle tight-binding parameters is the preliminary step in 
constructing more realistic many-body Hamiltonians. The role of the band theory is, 
thus, quite ambivalent: on one hand, it is the final "language" used in efforts towards 
understanding a broad variety of phenomena; on the other hand, it is the starting point 
in developing realistic interaction Hamiltonians for sophisticated phenomena such as 
magnetism and superconductivity. 

The tight-binding method is the simplest one employed in the electron band 
calculations and it is described in every textbook in solid state physics; the layered 
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perovskites are now probably the best investigated materials and the Fermi surface 
is a fundamental notion in the physics of metals. There is a consensus that the 
superconductivity of layered perovskites is related to electron processes in the Cu02 
and Ru02 planes of these materials. It is not, however, fair to criticise a given study, 
employing the tight-binding method as an interpolation scheme to the first principles 
calculations, for not thoroughly discussing the many-body effects. The criticism should 
rather be readdressed to the ab initio band calculations. An interpolation scheme cannot 
contain more information than the underlying theory. It is not erroneous if such a 
scheme works with an accuracy high enough to adequately describe both the theory and 
experiment. 

In view of the above, we find it very strange that there are no simple interpolation 
formulae for the Fermi surfaces available in the literature and the experimental data are 
being published without an attempt towards simple interpretation. One of the aims of 
the present paper is to help interpret the experimental data by the tight-binding method 
as well as setting up notions in the analysis of the ab initio calculations. 



2. Layered cuprates 

2.1. Model 

The Cu02 plane appears as a common structural detail for all layered cuprates. 
Therefore, in order to retain the generality of the considerations, the electronic properties 
of the bare Cu02 plane will be addressed without taking into account structural 
details such as dimpling, orthorhombic distortion, double planes, surrounding chains 
etc. For the square unit cell with lattice constant three-atomic basis is assumed 
{Rcu, Rob5 f^Obl = {O5 (<^o/25 0), (0, aQ/2)}. The unit cell is indexed by vector n = 
{ux^riy)^ where riy = integer. Within such an idealized model the LCAO wave 
function spanned over the |Cu3(ia;2_^2), |Cu45), \0s,2px), |0b2p^) states reads as 

n 

+^n^Cu4.(r - H-Cu - ^^0^) + i:^n^cu3d(r - Rcu " ^^o^)] , (2.1) 

where = (^n, X^^ Y^) is the tight-binding wave function in lattice representation. 

The neglect of the differential overlap leads to an LCAO Hamiltonian of the Cu02 
plane 

H = Y.{Dl[-tpd{-Xn + Xx-i,y + - Yx,y-i) + eaD^] 

n 

+ Sl^[~tsp{ — Xn + Xx-l^y — + Yx^y-i) + esSn] 

+ Xn[~'f^pp0^n ~ Yx^l^y — Yx^y-l + Yx^l^y-l) (2-2) 

~ tsp{~Sn + Sx-\-l,y) — tpd{ — Dn + Dx-\-l^y) + €pXn] 
+ ^l[~'tpp{Xn — Xx-l^y — Xx^y-\-l + Xx-l^y-\-l) 
~tsp{~Sn + Sx,y-\-l) — tpd{Dn + Dx^y-\-l) + 
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where the components of should be considered as being Fermi operators. The 
notations e^, e^, and stand respectively for the Cu3da,2_j^2, Cu4s and 02pa single- 
site energies. The direct O^^Px — ^ Oi,2py exchange is denoted by tpp and similarly tsp 
and tpd denote the Cu4s 02p and 02p Cu3da,2_y2 hoppings respectively. The sign 
rules for the hopping amplitudes are sketched in figure ^] — the bonding orbitals enter the 
Hamiltonian with a negative sign. The latter follows directly from the surface integral 
approximation for the transfer amplitudes, given in [Appendix A . 
For the Bloch states diagonalizing the Hamiltonian ( p.2|) 



Sn 



E 



e'^-^ATp 



^ipn 



(2.3) 



where is the number of the unit ceUs, we use the same phases as in references [jT2|, [13 

= ^{Px — 7t)^ Lph = ^{py — tt). This equation describes the Fourier transformation 
between the coordinate representation = (^n, S'n, -^n, ^n), with n being the ceU 
index, and the momentum representation ipp = {Dp^ Sp^ X^, Yp) of the TB wave function 
(when used as an index, the electron quasi-momentum vector is denoted by p). Hence, 
the Schrodinger equation ihdt^^^p^a = [4^p,ai H] for ipp^a{t) = ^~^^^^^i^p,a^ with a being the 
spin index (|, |) (suppressed hereafter), takes the form 



(//(^'^) - el) ijp 



where 



and 








^pd^ X 


tpd^Y ^ 


(Dp\ 







tsp^x 


tsp^Y 


Sp 


^pd^ X 


^sp^ X 


^p 


^PP^ X 


Xp 


tpd^Y 


tsp^Y 


tpp^X 




\Yp) 


- ^dl 




^S) ^p ^ 


^pi 





0, 



(2.4) 



= 2sin(|p^), = 2sin(|p^), x = sin^(|p^), y = sin^(|p^) 
< Px^ Py < 27r. 

This 4cr-band Hamiltonian is generic for the layered cuprates, cf. reference [|13[. We have 
also included the direct oxygen-oxygen exchange tpp dominated by the a amplitude. The 
secular equation 

det (//^^^) - eU) = Axy + B{x + y)+C = (2.5) 
gives the spectrum and the canonical form of the CEC with energy-dependent coefficients 



A{e) - 16(4t^^t^^ + 2tlptpp£d - 2tl^tpp£s - tlp£d£s) 



B{e) 
C{e) 



^^pitlpSd + tldSs) 



^d^s^p- 



Hence, the explicit CEC equation reads as 
Py 



Bx + C 

± arcsin y^, if < ?/ = — - — < 1. 

Ax + B 



(2.6) 



(2.7) 



Modelling of layered superconducting perovskites 



6 



This equation reproduces the rounded square-shaped FS, centered at the {tt^tt) point, 
inherent for aU layered cuprates. The best fit is achieved when B and C are considered 
as fitting parameters. Thus, for a CEC passing through the D = {pd'>Pd) and C = (pc, tt) 
reference points, as indicated in figure |2|, the fitting coefficients (distinguished by the 
subscript /) in the canonical equation Afxy + Bf{x + y) + Cf = have the form 

^/ = 2xd - - 1, Xd = sin^(f>d/2) 

Bf = x,-xl, x, = ^m\pj2) (2.8) 

Cf = xl{xc + 1) - 2XcXd, 

and the resulting LCAO Fermi contour is quite compatible with the LDA calculations 
for Nd2-a;Cea^Cu04_(5 [1^, |]. Due to the simple shape of the FS the curves just coincide. 
We note also that the canonical equation (|2.5|) would formally correspond to 1-band TB 



Hamiltonian of a 2D square lattice of the form 

e(p) = —2t{co8px + C08 Py) + cospx cos Py^ 

with strong energy dependence of the hopping parameters, where is the anti-bonding 
hopping between the sites along the diagonal, cf. references [|16|, W7 . 



2.2. Effective Hamiltonians 

Studies of the electronic structure of the layered cuprates have unambiguously proved 
the existence of a large hole pocket — a rounded square centred at the (tt, tt) point. This 
observation is indicative for a Fermi level located in a single band of dominant Cu3rfa;2_^2 
character. To address this band and the related wave functions it is therefore convenient 
an effective Cu-Hamiltonian to be derived by Lowdin downfolding of the oxygen orbitals. 
This is equivalent to expressing the oxygen amplitudes from the third and fourth rows 
of il 



X 

Y - 



o2 ^2 



(2.9) 



where r]p = Sp — -f^s^s'^^ and substituting back into the first and the second rows of the 
same equation. Such a downfolding procedure results in the following energy-dependent 
copper Hamiltonian 

+ + y + ^^^^ {x - y) \ 



Hcuie) 



Vp 



('2tpd)('2tsp) 



Vp 



(x-y) 



es+^^{x + y-'-^xy) ) 



(2.10) 



which enters the effective Schrodinger equation i/cu(^) = ^(^)- Thus, from (|2l9|) and 
( [j.l(J| ) one can easily obtain an approximate expression for the eigenvector corresponding 
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to a dominant CvL3dx2_y2 



character. Taking D 
the hopping amphtudes tw one has 



1, in the lowest order with respect to 



|Cu342_j^2) 



i.e. |Xp + |y|2 + |^|2 < 



S 
X 

\yj 



.) 



(2.11) 



/ 1 

{tpd/Vp)Sx 

-{tpdhp)sY 

1. We note that within this Cu scenario the Fermi 



V 



level location and the CEC shape are not sensitive to the tpp parameter. Therefore one 
can neglect the oxygen-oxygen hopping as was done, for example, by Andersen et al 
T2|, [l^ (the importance of the t^^ parameter has been considered by Markiewicz [|T4|1 ) 



and the band structure of the Hamiltonian (|2.10|) for the same set of energy parameters 
as used in reference [|13[ is shown in figure H (a). In this case the FS can be fitted by 
its diagonal alone, i.e. using only D as a reference point. Hence an equation for the 
Fermi energy follows, A{e^)x\ -\- 2B{e^)xd + C{e^) = 0, which yields = 2.5 eV. As 



IS 



seen in figure ^ (6), the deviation from the two-parametr fit, discussed in section 
almost vanishing thus justifying the neglect of tpp and using one-par ametr fit. 

However, despite the excellent agreement between the LDA calculations, the LCAO 
fit and the ARPES data regarding the FS shape, the theoretically calculated conduction 
band width Wc in the layered cuprates is overestimated by a factor of 2 or even 3 0. 
Such a discrepancy may well point to some alternative interpretations of the available 
experimental data. In the following section we shall consider the possibility for a Fermi 
level lying in an oxygen band. 



2.2.1. Oxygen scenario: the Abrikosov-Falkovsky model Various hints currently exist 
in favor of 02p character of the states near the Fermi level [|18|, [19|]- We consider that 
these arguments cannot be a priori ignored. This is best see if, following Abrikosov 
and Falkovsky [|20[ , the experimental data are interpreted within an alternative oxygen 
scenario. 

Accordingly, the oxygen 2p level is assumed to lie above the Cu3rfa^2_^2 level, and 
the Fermi level to fall into the upper oxygen band, Cd < Cp < < ^s- The Cvi3dx2_y2 
band is completely filled in the metallic phase and the holes are found to be in the 
approximately half-filled 02pa bands. To inspect such a possibility in detail we use 
again the Lowdin downfolding procedure now applied to Cu orbitals. From the first and 
second rows of ( p.4|) we express the copper amplitudes 

D = tif{s,X-s,Y) 

(2.12) 

S='if{s,X + s,Y) 

and substitute them in the third and the fourth rows. This leads to an effective oxygen 
Hamiltonian of the form 

i 



Ho{e) = B 



eff 



SyS^ 



(2.13) 
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with spectrum 



e{p) = 2B{e)ix + y) 



-1± Jl + (2r + r2 



(2.14) 



where 



B{e) 



+ i V 



^eff(e) = tpp + 2 



y-2 



T(e) = t^s/B 



(2.15) 
(2.16) 

corresponds to the sign for |r| < 1. 



and the conduction band dispersion rate ec(p) 
It should be noted that (|2.14| ) is an exact result within the adopted 4cr-band model. As a 
consequence, it is easily realised that along the (0, 0)-(7r, 0) direction the conduction band 
is dispersionless, ec{px^O) = 0. This corresponds to the extended Van Hove singularity 
observed in the ARPES experiment and we consider it being a hint in favour of the 
oxygen scenario (the copper model would give instead the usual Van Hove scenario). 

Depending on the r value two different limit cases occur. For r <C 1 one gets a 
simple Pade approximant 



6c (p) = 4teff(ec) 

and eigenvector of 
|c) = 



2xy 
X + y 



(2.17) 



S 
X 

\y J 



V 



(2.18) 



normalized according to the inequality + \S\'^ <C + |yp ^ 1. This limit case 
acceptably describes the experimental ARPES data e.g. for Nd2-xCexCn04-s^ material 
with single Cu02 planes and no other complicating structural details. Schematic 
representation of the energy surface defined by (|2.17| ) is shown in figure § (a). In figure || 
(6) we have presented a comparison between the ARPES data from reference and 
the Fermi contour calculated according to (|2.17| ) for x = 0.15. Note that no fitting 
parameters are used and this contour should be referred to as an ab initio calculation 
of the FS. 

The opposite limit case igff 3> i.e. r S> 1, has been analysed in detail 
by Abrikosov and Falkovsky [^. The conduction band dispersion rate 6c and the 
corresponding eigenvector of the Hamiltonian Hq (|2.13|) now take the form 

ec(p) = 4teff(ec)v^ (2.19) 

/ {tpd/£d){Sx +Sy)\ 



1 

72 



{tsp/£s){Sx 
1 
-1 



(2.20) 



/ 



provided that \D\^ + IS'p <C |Xp + |Fp ~ 1. In other words, the last approximation, 
r 3> 1, corresponds to a pure oxygen model where only hoppings between oxygen ions 
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are taken into account. Clearly, this model is the complementary to the copper scenario 
and is based on an effect completely neglected in its copper "counterpart" , where tpp = 0. 
This limit case of the oxygen scenario suitably describes the ARUPS experimental data 
for Pbo.42Bii.73Sri.94Cai.3Cui.9208+a: 0. The FS of the latter is fitted by its diagonal 
(the D point) according to the Abrikosov-Falkovsky relation ( |2.19|) and the result is 
shown in figure ^. 

There exist a tremendous number of ARPES/ ARUPS data for layered cuprates 
which makes the reviewing of all those spectra impossible. To illustrate our TB model 
we have chosen data for the Pb substitution for Bi in Bi2Sr2CaCu208, see figure |. 
In this case the Cu02 planes are quite flat and the ARPES data are not distorted by 
structural details. When present, distortions were misinterpreted as a manifestation 
of strong antiferromagnetic correlations. We believe, however, that the experiment by 
Aebi et al [|| reveals the main feature of the Cu02 plane band structure — the large 
hole pocket found to be in agreement with the one-particle band calculations. 

Besides the good agreement between the theory and the experiment, regarding 
the FS shape, we should also point out the compatibility between the calculated 
and the experimental conduction bandwidth. Indeed, within the Abrikosov-Falkovsky 
model pO[, according to ( |2.19|) , one gets for the conduction bandwidth < ^c(p) ^ '^c ^ 



Atpp^ which coincides with the value obtained from ( |2.17|) provided that t^^ <C tpp{e^—ed)- 



The ab initio calculation of tpp as a surface integral (see [Appendix A|) , making use 



of atomic wave functions standard for the quantum mechanical calculations, gives 
tpp ^ 200-350 meV in different estimations. This range is in acceptable agreement 
with the experimental t(;c — 1 eV [Q; within the LCAO model an exact analytic result 
for Wc can be obtained from the equation Wc = ^tpp + 8t^^/ {w^^ — Cd). 

We note also that the TB analysis allows the bands to be unambiguously classified 
with respect to the atomic levels from which they arise. Within such terms, for the 
oxygen scenario one can describe the metal^insulator transition as being the charge 
transfer Cu^+02^ C\i^^Ol~. The possibility for monovalent copper Cu^+ in the 



superconducting state is discussed, for example, by Romberg et al p2 . 
3. Conduction bands of Ru02 plane 

Sr2Ru04 is the first coper-free perovskite superconductor isostructural to the high- 
Tc cuprates The layered ruthenates, just like the layered cuprates, are strongly 



anisotropic and in a first approximation the nature of the conduction band(s) can be 
understood by analysing the bare Ru02 plane. One should repeat the same steps as 
in the previous section but now having Ru instead of Cu and the Fermi level located 
in the metallic bands of RnAdn character. To be specific, the conduction bands arise 
form the hybridisation between the RnAdxy^ Ru4rf^^, RnAdzx and Oa2p^, 0]j2px^ Oa,b2pz 
TT-orbitals. The LCAO wave function spanned over the four perpendicular to the Ru02 
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plane orbit als reads as 

+ e'^'^Za^^ij^^,^^ (r - Ro, - a„n) + e''^'Z6,„^„^,^ 



(3.1) 



r - Roi, - 



^ip n 



Hence, 


the TT-analog 


of (2.4) takes the form 














tz,zx^x 


0\ 








- el) tAj^) = 




'^Z,ZX^ X 
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Dzy 
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tz^zy^Y tzz^x^Y 


—£zb J 
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^ ^zxi 


^za ^ ^za? 


= 2cos(pa;/2), 








^ ^zy) 


•^2:6 ^ ^zb) ^Y 


= 2cos{py/2), 





0, 



(3.2) 



(3.3) 



and Czx^ ^zyi ^zai and e^b are the single site energies respectively for RnMzx^ Ru4rf^^ and 
0^2pz^ 0]j2pz orbitals. tzz stands for the hopping between the latter two orbitals and, if 
a negligible orthorhombic distortion is assumed, the metal-oxygen vr-hopping parameters 
are equal, tz^zy = tz.zx and also = = ^zb- The phase factors e^"^"'^ in are chosen 



in compliance with reference [|T3|], see equation ( U^ - 

Identically, writing the LCAO wave function spanned over the three in-plane tt- 
orbitals RnAdxy^ Oa2]9^, and 0\^2px in the way in which is designed one has for the 
"in-plane" Schrodinger equation 



f -8 
' ^xy 

tpdiT^x 

Y tpdlT^Y 



tpdn^x 
~^ya 



tpdiT^^ 

^pp^X^Y 



\ 



Ya 



= 0, 



(3.4) 



where tpdir denotes the hopping RuAdxy Oa,b2p7r and t^^, respectively,© a 2]9^ 0\^2px- 
The definitions for the other energy parameters are in analogy to ( p.3|) (for negligible 
orthorhombic distortion Cya = e^b 7^ ^z)- Thus, the vr-Hamiltonian of the Ru02 plane 
takes the form 



(3.5) 



p,a=T,i 



In a previous paper |23[ we have derived the corresponding secular equations and now 
we shaU only provide the final expressions in terms of the notations used here 

det(//("'^^) - el) = A^^'^'y^xy + ^(^'^^^(x + ?/) + C^"'^^) = 0, 



z,zx 



^zz^zx 



^^^zz^zx 



Z.ZX^ZX^Z) 



^zxi^z ^^^zz) 



J[{xy) _ S2tpptp^^ 

jg{xy) 
Q(xy) 



1 f\;^ 

^^^xy^pp 



— f^ F 
^pdiT^ya 

2 

^xy^ya 



(3.6) 



The three sheets of the Fermi surface in Sr2Ru04 fitted to the ARPES data given by 
Lu et al are shown in figure ^ (6). To determine the Hamiltonian parameters we 
have made use of the dispersion rate values at the high-symmetry points of the Brillouin 
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zone. To the best of our knowledge, the TB analysis of the Sr2Ru04 band structure 
was first performed in reference [|23[ (subsequently, the latter results were reproduced 
in reference [^] without referring to reference [|23[). The Ru02-plane band structure 



resulting from the set of parameters 

tzz = t^p = 0.3 eV, 8, = -2.3 eV, e,y = -1.62 eV, 

(3.7) 

tpdiT 1 eV, = -1.3 eV, Sya^^b -2.62 eV. 

is shown in figure ^ (a). This fit is subjected to the requirement of providing as good as 
possible a description of the narrow energy interval around whereas the filled bands far 
below the Fermi level match only qualitatively to the LDA calculations by Oguchi 0] and 
Singh [|] . In addition we note that the de Haas- van Alphen (dHvA) measurements 
of the Sr2Ru04 FS diflFer from the ARPES results 0. Thus, fitting the dHvA data 
by using modified TB parameters is a natural refinement of the proposed model. We 
note that the diamond-shaped hole pocket, centred at the X point (see figure ^ (&)), is 
very sensitive to the 'game of parameters'. For that band the van Hove energy is fairly 
close to the Fermi energy. As a result, a minor change in the parameters could drive a 
van Hove transition transforming this hole pocket to an electron one, centred at the F 
point. Indeed, such a band configuration has been recently observed also in the ARPES 



revision of the Sr2Ru04 Fermi surface |2J]. This can be easily traced already from the 
energy surfaces e{p) calculated earlier in reference [^]. The comparison of the ARPES 
data with TB energy surfaces could be a subject of a separate study. 

4. Discussion 

The LCAO analysis of the layered perovskites band structure, performed in the 
preceding sections, manifests a good compatibility with the experimental data and the 
band calculations as well. Due to the strong anisotropy of these materials, their FS 
within a reasonable approximation is determined by the properties of the bare Cu02 or 
Ru02 planes. 

Despite these planes having identical crystal structure, their electronic structures 
are quite different. While for the Ru02 plane the Fermi level crosses metallic 7r-bands, 
the conduction band of the Cu02 plane is described by a a-Hamiltonian ([2^ ). The 
latter gives for the Cu02 plane a large hole pocket centered at the (7r,7r) point. Its 
shape, if no additional sheets exist, is well described by the exact analytic results 
within the LCAO model. Equation ([2?5| ), as found for Nd2-xCea^Cu04_(5 [||, ^ and 
Pbo.42Bii.73Sri.94Cai.3Cui.9208+x p. For a number of other cuprates, YBa2Cu307_(5 p7|] . 



YBa2Cu408 0, Bi2Sr2CaCu208 |2|, Bi2Sr2Cu06 H, the infinite-layered 



superconductor Sri_a;Caa;Cu02 HgBa2Ca2Cu308+(5 [^], HgBa2Cu04+(5 

HgBa2Ca,_iCu,02.+2+c5 [0^ Tl2Ba2Ca,_iCu,04+2n, 0, Sr2Cu02F2, Sr2Cu02Cl2, 
Ca2Cu02Cl2 this large hole pocket is easily identified. For all of the above 

compounds, however, its shape is usually deformed due to appearance of additional 
sheets of the Fermi surface originating from accessoires of the crystal structure. 
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As the most important implication for the Cu02 plane we should point out 
the intrinsic alternative about the Fermi level location (see section g). It is 
commonly believed that the states at the FS are of dominant Cu3rfa;2_^2 character (see 
e.g. reference [|13|). Nevertheless, the spectroscopic data for the FS can be equally well 
interpreted within the oxygen scenario, according to which the FS states are of dominant 
02pa character. A number of indications exist in favour of the oxygen model and the 
importance of the tpp hopping amplitude 0, [l^ : 



(i) 01s 02p transitions observed in EELS experiments for the metallic phase of the 
layered cuprates, which reveal an unfilled 02p atomic shell; 

(a) the oxygen scenario reproduces in a natural way the extended van Hove singularity 
observed in the ARPES experiments while the Cu scenario fails to describe it; 

(Hi) the metal-insulator transition can be easily described; 

(iv) the width of the conduction band is directly related to the atomic wave functions. 



Some authors even "wager that the oxygen model will win" [|T9| (if the oxygen scenario 
is corroborated, due to the cancellation of the largest amplitude tgp the small hoppings 
tpd and tpp should be properly evaluated eventually as surface integrals (see |Appendix| 
^ and some band calculations may well need a revision) . It would be quite valuable if 
a muffin-tin calculation of H2" ion was performed and compared with the exact results 
when the hopping integral is comparatively small, of the order of the one that fits 
the ARPES data tpp ^ 200 meV. We also note that even the copper model gives an 
estimation for tpp closer to the experiment than the LDA calculations. The smallness 
of tpp within the oxygen scenario, on the other hand, is guaranteed by the nonbonding 
character of the conduction band. This scenario, therefore, can easily display heavy 
fermion behaviour, i.e. effective mass meff huge, and density of states (DOS) 

oc meff oc I /tpp (we note that no realistic band calculations for heavy fermion systems 
can be performed without employing the asymptotic methods from the atomic physics). 
It is also instructive to compare the TB analyses of heavy fermion systems and layered 
cuprates. The alternative for the Fermi level location (metallic vs. oxygen band) exist 



for the cubic bismuthates as well [^, When the Fermi level falls into heavy 



fermion oxygen bands, one of the isoenergy surfaces is a rounded cube [^]. Indeed, 
such a isoenergy surface has been recently confirmed by the LMTO method applied to 
Bao.6Ko.4Bi03 H- 



Due to the equally good fit of the results for the FS of the layered cuprates within 
the two models we can infer that at present any final judgement about this alternative 
would be premature. Thus far we consider that the oxygen model should be taken 
into account in the interpretation of the experimental data. Moreover, the angular 
dependence of the superconducting order parameter A(p) oc co8{px) — co8{py) is readily 
derived within the standard BCS treatment of the oxygen-oxygen superexchange |36|] . 
Analysis of some extra spectroscopic data by means of different models would finally 
solve this dilemma. This cannot be done within the framework of the TB method. A 
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coherent picture requires a thorough study, where the TB model is just an useful tool 
to test the properties of a given solution. 

Up to now, the applicability of the LCAO approximation to the electron structure of 
the layered cuprates can be considered as being proved. The basis function of the LCAO 
Hamiltonian can be included in a realistic one-electron part of the lattice Hamiltonians 
for the layered perovskites. This is an indispensable step preceding the inclusion of 
the electron-electron superexchange, elect ron-phonon interaction or any other kind of 
interaction between conducting electrons. 
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Appendix A. Calculation of O-O hopping amplitude by the surface integral 
method 



From quantum mechanics ^Tj it is well known that the usual for the quantum chemistry 
calculation of the hopping integrals as matrix elements of the single particle Hamiltonian 
does not work when the overlap between the atomic functions is too weak. If the hopping 
integrals are much smaller than the detachment energy, they should be calculated as 
surface integrals using (eventually distorted by the polarisation) atomic wave functions. 

Such an approach has been applied by Landau and Lifshitz and Herring 
and Flicker [|38|] to the simple H2" problem and now the asymptotic methods are well 
developed in the physics of atomic collisions [^. On the basis of the above problem 
one can easily verify that the atomic sphere muffin-tin approximation of the Coulomb 
potentials usual for the condensed matter physics undergoes fiasco when the hopping 
integrals are of the order of 200-300 meV. Therefore, the factor 2-3 misfit for a single 
electron problem cannot be ascribed to the strong-correlation effects, renormalizations 
and other incantations which are often used to account for the discrepancy between the 
experimental bandwidth and the LDA calculations. 

Usually condensed matter physics does not need asymptotically accurate methods 
for calculation of hopping integrals which leads to zero overlap between the muffin-tin 
and asymptotic methods. However, for the perovskites the largest hopping tgp cancels in 
the expression for the upper oxygen band ec(p). Thus, small hoppings become essential, 
but having no infiuence on the other bands, and the necessity of taking into account the 
tpp is of topological nature. 
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Following the calculations for H2" in a simplified picture of two oxygen atoms 
Oa, Ob separated by distance d = the surface integral method gives for the oxygen- 
oxygen exchange the following explicit expression 

tpp = ^ II ^^oA^o^ - ^o^^oM^^y^ (A.i) 

s 

where the integral is taken over the surface S halving and m is the electron mass. Thus 
tpp = tpp{^) is a function oi ^ = /^|Roa ~ f^Obl with /^^/2 being the oxygen detachment 
energy in atomic units and the detailed derivation of ( |A.1|) can be found, for example, 



in reference [39 



We note that the derivation of tpp{^) imposes no restrictions on the basis set 
{ip} used. Hence we choose {ip^ ^} to be the simplest minimal (MINI) basis 
used EOl, for example, in the GAMESS package for doing ab initio electronic structure 



calculations |4l|]. The MINI bases are three Gaussian expansions of each atomic orbital. 
The exponents and contraction coefficients are optimised for each element, and the s 
and p exponents are not constrained to be equal. 

Accordingly, the oxygen 2p radial wave function R2p{r) is replaced by a Gaussian 
expansion R2^\r) and has the form 

4?(C,r) = j2C2p,i92pAC2p,i,r), (A.2) 



where g2p{Ci ^) = ^2p,i^~^^^''^ , and the coefficients for oxygen are given in table |T]. It is 
then normalized to unity according to R^^'^^r'^dr = 1. 

By multiplying with the corresponding cubic harmonic the oxygen wave functions 
are brought into the form 



,«„.(rj = <'(C, r.),/i-i!i, ( = ; ^ , (A.3) 

^ \/ 47r Ta Ta = rJ 



and analogically for ip^ (r — Rob)- Substituting (|A.3|) in ( |A.1| ) we get 



C™' = 340 meV. 



In reference the same integral has been calculated with {ip} being the asymptotic 
wave functions ||S^ appropriately tailored to the MINI basis at their outermost inffection 
points r^^\ i.e. 

R,p{r) = (A.4) 

\ r ' — 

with tv = 0.329 and A = 0.5. The value obtained is 

^(asymp) _ 2IO mcV, 

found to be in good agreement with that fitted from the ARPES experiment within the 
oxygen scenario. Similar calculation, for example, gives for the tpd and tgp hoppings 

(MINI) (MINI) 

tp^ = 580 meV, t^^ ^2.5 eV. 
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Note added in proof. In a very recent paper by Campuzano J C et al 1999 [Phys. Rev. Lett. 
83 3709] the ARPES Fermi surface of pure Bi2Sr2CaCu208+5 has been presented in the inset 
of their figure 1 (a). This experimental finding is in excellent agreement with our tight-binding 
fit to the Fermi surface of Pbo.42Bi1.73Sr1.94Ca1.3Cu1. 9208+0:5 studied by Schwaller P and co- 
workers in reference 6, given in figure 5 of the present paper. The remarkable coincidence of 
the Fermi surfaces of these two compounds is a nice confirmation that Pb substitution for Bi 
is irrelevant for the band structure of the Cu02 plane and the Fermi surface of the latter is 
therefore revealed to be a common feature. 
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Figure 1. Schematic of a Cu02 plane (only orbitals relevant to the discussion are 
depicted). The solid square represents the unit cell with respect to which the positions 
of the other cells are determined. The indices of the wave function amplitudes involved 
in the LCAO Hamiltonian ( |2.2|) are given in brackets. The rules for determining the 
signs of hopping integrals tpd^ tsp, and tpp are shown as well. 



Figure 2. LDA Fermi contour of Nd2-a?Cea;Cu04_5 (dotted line) calculated by Yu 
and Freeman jij (with the kind permission of the authors), and the LCAO fit (solid 
line) according to (|2.5| ). The fitting procedure uses C and D as reference points. 



Figure 3. (a) Electron band structure of the generic for the Cu02 plane 4<j-band 
Hamiltonian using the parameters from reference [|l^ and the Fermi level = 2.5 eV 
fitted from the LDA calculation by Yu and Freeman 0; (&) The LCAO Fermi contour 
(solid line) fitted to the LDA Fermi surface (dashed line) for Nd2-xCexC\i04-s 
using only D as a reference point. The deviation of the fit at the C point is negligible. 



Figure 4. (a) Energy dispersion of the nonbonding oxygen band ec(p), 

equation ( ^.17| ). A few cuts through the energy surface, i.e. CEC, are presented 
together with the dispersion along the high symmetry lines in the Brilloiun zone; (b) 
The Fermi surface of Nd2-a^Cea;Cu04_<5 (solid line) determined by equation ( 2.17 ) for 
X = 0.15 (shaded slice in panel (a)) and compared with experimental data (points 
with error bars) for the same value of x after King et al and cp denote the polar 
and azimuthal emission angles, respectively, measured in degrees. The empty dashed 
circles show k-space locations where ARPES experiments have been performed (cf. 
figure 2 in reference 0) and their diameter corresponds to 2° experimental resolution. 



Figure 5. (a) ARUPS Fermi surface of Pbo.42Bii.73Sri.94Cai.3Cui.9208+x by Aebi 
et al (b) LCAO fit to (a) according to the Abrikosov-Falkovsky model using 
the D reference point with = 0.171 x 27r. 



Figure 6. (a) LCAO band structure of Sr2Ru04 according to (|3^). The Fermi level 
(dashed line) crosses the three KuAde bands of the Ru02 plane; {b) LCAO fit (solid 
lines) to the ARPES data (circles) by Lu et al [Q, cf. also reference | p3[ . 
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Table 1. Coefficients for tiie oxygen 2p wave function in ttie MINI basis |4l|] . 



i 








1 


8.2741400 


2.485782 


0.708520 


2 


1.1715463 


1.333720 


0.476594 


3 


0.3030130 


0.263299 


0.130440 
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